HAS Iffacac)* Icport Mo. 146 


May 1973 


THEBMAL EMGDIEERIHG RESEARCH 





Csci 


^Of; 







Thla research work was supported 
The KASA/MSFC 

Under UAH-MSFC/NASA Cooperative Agreeaent 


by 

Modification Mo. 


1 


The University of Alabaaa in Huntsville 
Huntsville, Alabaaa 



UAH Rese«ch Report No. 146 


May 1973 


THERMAL ENGINEERING RESEARCH 


Final Report 


This research work was supported by 
The NASA/MSFC 

Under UAH-MSFC/NASA Cooperative Agreement Modification No. 1 


The University of Alabama in Huntsville 
Huntsville, Alabama 



PREPACK 


This final report of efforts for the UAH-MSFC/NASA Cooperative Agreement 

Modification No. 1 (Thermal Engineering Research) presents the results of the 

following projects for the period from September 1971 through May 1973: 

PROJECT 1 Analytical Treatment of Gas Flows Through Multilayer 
Insulation 

PROJECT II Study of Deposition Model of Frost and Ice From Humid Gas 
in Ducts 

PROJECT III Analysis of Thermal Environment in the Thrustor Cavity of 
Space Vehicles 

Due to the diversity of the projects* requirements, the projects were assigned 
to three teams of research personnel. Under the general supervision of the 
Principal Investigator, the responsibility of technical direction for Project 
I was assigned to J. T. Lin; Project II was originally assigned to R. W. 
Blanton, and subsequently to C. C. Shih; Project III was assigned to J. J. 
Brainerd. Project I was carried out solely by J. T. Lin with minor assistance 
from the F&TE Lab staff. The bulk of work for Project II was carried out by 
B. K. Singh, then Bhuminder Singh. Project III was conducted mainly by Jim 
Chiou under proper technical guidance. 

To make the success of this study possible valuable assistance was pro- 
vided by the following staff members of the Fluid Dynamics Laboratory at the 
F&TE Laboratory of The University of Alabama in Huntsville: 

Margaret P. Cuthlll Carol Holladay Debbie Maples 

Tommie McMeans Don Miller Barbra Oedamer 

Thomas A. White and other undergraduate lab assistants. 

The computer laboratory, instrumentation laboratory, and machine shop of 
The University of Alabama in Huntsville also provided useful assistance to 
the projects. 


Cornelius C. Shih 
Principal Investigator 
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ABSTRACT 


A thcoreHcal investigation of gas flow inside a multilayer insulation system 
has been mode for the case of the broadside pumping process. A set of simultaneous 
first-order differential equations for the temperature and pressure of the gas mixture 
was obtained by considering the diffusion mechanism of the gas molecules through 
the perforations on the insulation layers. A modified Runge-Kutta method was used 
for numerical experiment. The numerical stability problem was investigated. It 
has been shown that when the relaxation time is small compared with the time period 
over which the gas properties change oppreciably/ the set of differential equations 
can be replaced by a set of algebraic equations for solution. Numerical examples 
were given^and comparisons with experimental data were mode. 
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Subscripts: 


Superscripts: 


area of one insulation layer (one side) 

total area of perforation on one insulation layer 

a' = (1 - c) A 

perforation coefficient c ~ a/A 
Boltzmann constant 

number of interstitial spaces between insulation layers 

number of kinds of absorbed gases plus 1 

mass of a gas molecule or an integer 

number of gas molecules 

number density of gas molecules or an integer 

pressure 

outgassing rates, see Equations (4), (7), and (12) 

temperature 

time 

volume 

mean thermal velocity of gas molecules 
number of degrees of freedom of a gas molecule 
.elaxation time 

subscript c refers to quantities inside the vacuum chamber; f 
to condition of the fuel tank; a to atmospheric condition; other 
subscripts refer to quantities in a certain interstitial space 
superscript 1 refers to the purge gos, other superscripts refer 
to a certain kind of gas molecules 



1 . Introduction 


Recent developments in propulsion technology hove stimulated interest in 
the study of multiloyer iniulation systems for the fuel tank of a rocket booster. The 
purpose of multilayer insulations is to eliminate beat conduction between the fuel 
tank and the environment. It is then desirable that the spaces between the insula- 
tion layers be mointoined at low pressures. However, it is not practical to require 
complete vacuum in the spaces between the insviation layers, although it would be 
ideal in principle, for the reasons that^ complete vacuum is hardly attainable^and 
that when the pressure gets too low, deformation of the layers would cause solid 

I 

contacts/ which in turn would induce hect conduction. It is therefore customary 

-4 

in practice to require a pressure to be of the order of 10 torr. This low pressure 

is usually achieved before launch from the ambient pressure by some pumping device. 

There hove been two pumping procedures in practice; namely, the broadside pumping 

by which the direction of gas flow is perpendicular to the insulation layers, and 

the edge pumping which cousesthe gas to flow parallel to the layers. Experiments 
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of both evacuation procedures have been performed and results reported. ' ' The 
present work is a theoretical analysis for the broadside pumping process. 

In the case of broadside pumping, the insulation layers are perforated so 
that gas molecules can go through small holes on the layers, resulting in streaming 
gas motion perpendiculor to tne insulation layers. It is obvious that, in order to 
minimize layer deformation and solid contacts, small pressure differentials across 
the layers should be maintained during evacuation. This con be achieved by small 
perforation, i.e., the total area of the holes on an insulation layer is small com- 
pared to the total area of that layer. If we further require that the diameter of a 
single holebesmall or comparable to the distance between two successive layers, 
there is the advantage of uniformity of gas motion, which renders convenience dur- 
ing operation. 

A multilayer insulation system consists of a large number^ 0^ ^he order of 10 ) of 
extremely thin sheets of low thermal conductivity and lightweight materials, kept 
parallel to one another with a total thickness of about 1". 



2* Formulation 


Consider a fuel tank of simple geometry wrapped with multilayer insulation 
and situated inside a vacuum chamber. Consider the broadside pumping process 
and assume that the perforation areas are small . Then, a gas molecule in any inter- 
stitial space between two successive layers will# on the average# collide many times 
with the walls or with some other gas molecules in that interstice before getting 
through a hole to a neighboring interstice. Thus, it is plaussible to assume that the 
gas in any interstitial space is in thermal equilibrium with the temperature and 
pressure pertinent to that interstice# and the motion of the gas is simply a diffusion 
process* 

Consider the gas inside the insulation system as a mixture of a certain kind 
of purge gas and a number of different kinds of gases originally absorbed in the in- 
sulotion materials. Let be the total area of the outermost insulation layer (one 
side)/ and c^ the perforation coefficient of the same layer (the ratio of the total 
perforation area on the first layer to A.^)# and so on. Let Vj be the volume of 
the interstice between the first and the second insulation layers# and Nj the number 
of molecules in the volume of the Interstice between the second and the 

third insulotion layers/ and N 2 the number of molecules in y and so on. Since 
the gases in any Interstice are in thermal equilibrium# each component gas has a 
Maxwellian distribution with its local density ond the common local temperature. 

From the kinetic theory of gases# we have# for the gases in the i th interstice# 
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where o, “ ®j ^ °'j ~ ^ ^ 1^® Boltzmann constant, 6 the num- 

ber of degrees of freedom of motion of a gas molecule, n., v. and T, are the number 
density, the mean thermal velocity and the temperature of the gas molecules in the 
i th interstice, respectively. The superscript j refers to the purge gas if j ~ 1/ 
or the absorbed gas of kind j otherwise. The number of interstitial spaces is L, 
ond the number of kinds of the absorbed gases is (M-1) so that 
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Finally, 



Ecjuation (1 )hi I'lti's the 


0 (j =1) 

(4) 

number of outgassing molecules/ of kind j from the 
walls binding the i th Interstice, per unit area and 
time ( j ^ 1 ). 

conservation of nuiss and Eqs. (2) the conservation of energ 
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since the L th interstice IS n^.*xL to the fuel tank, we may assume that 

0 ^ ^ I * Oand temperature of the fuel tank. Since inter- 

stice 1 is the outermost/ denoting the number densities, the mean thermal velocities 
end the temperature of the molecules in the vacuum chamber by v ^^and 

T / respectively/ we have n ^^ = n v v ^^andT ~T . Notice that 

' O C O C...OC 

as the volume of the vacuum chamber is very large, n (j 1) are increasingly 

c 

smaller compared to at all times and wc may assume that n^ ^ ^ ^=0(jf^l) . 

By making this assumption/ v (j 1 )btM omes meanlnRless and it drops out Lrom 

Eqs. (l)and (2) automatically . Therefore/ in Eqs. (l)and (2), we use 
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0 , T = T , T «T , V V 
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n 0) ^ (i)«Q 


Following p. ~ ^ T. / where p. is the pressure in the ith interstice/ 


we moy assume that 
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p.= E p.fi> 

III I TTr^ M 


( 6 ) 


which is plaussible when the pressure is not too high. 

Using Eqs. (3) and (6), il tl»e out^assin^ rat(*s arc expressed in the 

customary units of pressure times volume per unit area and time denoted by 


= kT. 

I I 


(7) 


then/ Eqs. (1)/ (2)/ and (5) can be written as 
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( 10 ) 


p p and p = 0 (j 7^1) . 

c ’^o 
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where p ^ ^ ts the chamber pressure due to the purge gas. 
c 

The mean thermal veiccities are given by 


/ g L T \ 


-,M. (11) 


where m^^ is the mass of a molecule of the j th kind. 


(j) 


To obtain a solution to the problem/ the outgassing rates Q. must be 

1 * 

known a priori / usually determined experimentally. ' * Outgassing is a process 

In which gas molecules originally adsorbed or absorbed In a solid material leave 

that material under reduced pressure or elevated temperature. The outgassing rate 

of the insulation layers depends on the material and preconditioning of the layers, 

the temperature, Che pri^f^suro, the time nnd the kinds ot absorbed : a 

specific material and preconditioning and a sporific absorbed gas, Q; ' 

(T,p,t). Around the room temperat tire ran^e, tlier<‘ is not outgassing at atmospheric 
or hi gher pressures . Outgass i ng oc eur.s when t he bnrkp”onnd pressure is reduced 
considerably below tht* atmosplierlc pressure. Tlu* outgassing rate increases as the 
pressiire decreases and attains apprcclblt‘ values only at very low pressures. 
However, In very low |)ressure ranges, the outgassing rate is a weak function of 
pressure. Therefore, Its dependence on the pressure may be approximated by 
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where p. ~ .■4K / P is the atmospheric pressure and is theout- 

I J - 1 I a 

gauing rates at extremely low pressures. 


*** 


The absorbed gases are mainly water vapor, with small amounts of COj and 

Nj (Ref. 7). 
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The initial conditions associated with the system of Equations (8^ ond (9) 
are prescribed according to a specific situation encountered. Usually the evacuation 
process starts with atmospheric pressure inside the insulation system. In this case, 
the initial corditions can be written os 

Pj^^^ (o)= p^, P.^^o)= 0(j?^l),i = 1, 2, ---L 

(13) 

T.(0)=T 0/ i = 1/ 2, , L - 1. 

Now, with the (^ 1)L ~ 1 initial conditions given" by Eqs. (13) or pre- 
scribed otherwise,and wi th the om K rates known, and by using Eqs. (10) 

and (11), llien Eqs. (8) and (9) <-onstllntt‘ an apfiropr iat i- sot ;>f (Mi-1 ) 1,-lsimul- 
tanoouR first order dlfforentlal equations for tiu- same number of unkr.own 

functions of time, namely, Tj^(t)k -L, 2, — b-1 and n^^ i ^ (t) , i=l 1.; j=l, — ,M. 

In some engineering applications, e.g., the Apollo Telescope Mount insula- 
tion, the temperature is quite uniform across the entire insulation system, and the 
gas flow inside the system is nearly isothermal . In this case, the problem reduces 
to the solution of Eqs. (8) with Eqs. (10) and (II), and the initial conditions (13), 
when all temperatures are set equal to a constant. 
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3. Numerical Analysis 


In a nu'lHloyer insulation system, the number L usually is of the order of 

2 

10 • Thus, the governing Eqs* (8) and (9) consist of a large number of first order 
differential equations which ore nonHineor and coupled. The analytical solution 

of this system of equotions can hordly be obtained^ ilieri toro, om* a dirt-ct 

numi^rical mothod for its solution. 


For simplicity in numerical experiments, we took on isothermal case, and 

assumed water vapor was the only obsorbed gas present ( M = 2 ) • Furthermore, 

( 2 ) 

a hypothetical R' (T , t)and chamber conditions were used. When the tempera hjre 

o 

is constant, the energy equotions (Eqs* (9) )are not needed, and the problem reduces 
to the solution of Eqs* (8) with initio! conditions (13). In the process of numerico! 
computation, however, the problem of instability occurred. A numerical program 

g 

using a modified Rurige^Kutta method has been set up for this Initial value problem. 

Computer experiments on this scheme showed that the system was stable when the 

values of = (4 V. /a. ) (^ m / 8 k T ) , which have a dimension of 

III o 

the time, were large* However, in order to establish numerical stability for small 
/• \ 

values of , the step sizes of integration had to be so small that a solution could 

not be obtained with o reosor)oble length of computer time* Since in ordinary ap* 

/. \ 

plications the values of are very small, we need to seek some olterrKjtive or 

approximation to the system of Eqs* (8), (9) and (13). 

For the case j ~ I, Eqs. (8) can be written as 
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time. 

i - 1 
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1 /T, I 

1 + I 1 I 1 
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dimensionless quant! Hes 


or order 1 . 
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Let the lorgest of oli the four defined above be t • Define the non~ 


m 


dimensionoi temperatures or>d pressures os T'. ~ T. /T,, where T iso reference 

' ' 0) ^ 

temperature the room temperature), ond p'j ~P j • 

period over which T'. ond p! chortge appreciably, and define the non-dimensional 
time os t’-t/r so thot dT*. /dt'ond dp'. /dt' hove ULiKniiiiJos *>;■ 4>rii. 7 

C I ^ t 

1- 

In lerms of T*.# p\ and t* , Eqs. (14) can be wrirten as 
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where $. • T / T. efc-, ore <fi mens! onl ess quantities of order 1. 

I m I 

We sec that 7 is a measure of the relaxation time- If 7 < < t ^ fhe 

m me 

two terms on the left hand side of Eqs- (15) ore vanishingly small compared to any 

term on the right hand side/ since the mognitudes of tK; two terms inside the brackets 

on the left bond side are of order 1 . In this case, component 1 o^ the gas mixture 

in the i th interstice is ’"quosi^teody”'/ which means that at ony insfont component 

1 in interstice i con be considered instantaneously steady and the derivative terms 

in its go ^eming equation can be dropped- This argument applies to al! other com~ 

ponents of the mixture in every interstice, i-e-, .ipplifs cquailv 1 1 to each 

equation of Eqs- (8) for coses j ^I and Eqs- (9) individually* 

Asa finol step towards completion of this argument, let 1 ^ m ^ L ond 

^n 1 f* ) 

1 < n ^ M be two integers such that 7 ' ' is the largest of oil " , ' defined by 
— m 1 
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Then, the reloxotion time of the gos mixture in the entire insulation system t 

is bounded by L t , ■ .e., 
vn 


T < (17) 

S ^ Tti 

Denoting the time period over which the properties of the gas mixture chonge 
appreciably by we can drop all the derivative terms in Eqs- (8) and (9) 
simultoneously if 

T or L T « T . (18) 

s in I 

Thus, this system of differential equations reduces to o set of algebraic 
equations* 

~2 

Now, for ordirwry engineering applications, ' is of the order 10 ~ 

“1 . • * 

10 sec., while is measured by hours. Therefore, condition (18) is usually 

fulfilled, and Eqs. (8)arKl (9) reduce to 



1 - 1, 2, , L - 1. 

io 


( 20 ) 



As Eqs. (19) and (20) ore cigebinic equations/ there is no need for the pre- 
scription of initial conditions* 

To obtain the solution for o particular problem, we first choose a sequence 

of the time (t., t^, — t )• The solution at any chosen time, soy tw is ob- 
I z n I 

tained by solving Eqs. (19) and (20) with the chamber conditions and fuel tank 
temperature at t^ (see Eqs. (10) )and \[ , t^ (see Eq. (12) ) which are 
pertinent to that problem. After we finish with the time sequence, we have the 
pressures and temperatures in oli the interstices of the insulation system os functions 
of the time. If the process is isothermal, one only has to solve Eqs. (19) with the 
relevant chamber conditions and outgassing rates. 
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4. Numerical Results and Comparison with Experimental Data 

For the purpose of comparison with experiment, we choose examples,the 
experimental data of which are available* 

A circular disk of 1“ thick and 6“ in diameter is composed of a variety of 
numbers of insolation sheets of crinkled single-alominized mylar. The edge of the 
disk is sealed with a solid insulation sheet (i .e., without perforation). The last 
loyer (the back sheet) is also solid, but the rest of the insulation layers are per- 
forated with a perforation coefficient equal to 0.0138. The perforation holes are 
0.09375” in diameter. The distance between two adjacent holes is 0.707" center 
to center . 

The disk is placed in a large vacuum chamber. The chamber and hence 

the multilayer system 0 *^be disk) are maintained at room temperature all the 

time, so it is on isothermal process. At the beginning of the experiment, the 

vacuum chamber and the multilayer system are filled with nitrogen ^urge gas), and 

the pres*'.’ e inside each interstitial space between two successive insulation layers 

end the pressure in the chamber are atmospheric. The chamber is then evocuoted. 

The pressure in the chamber and the pressure in the last interstice of the insulation 

system are recorded as time proceeds. 

/• \ 

The outgassing rotes, R^' as functions of the time, of the insulation sheets 

wore ob.ni.ied before hand by separate experiments. The outgas components are 

moir'* water vapor, with small amounts of nitrogen and carbon dioxide. For sim- 

( 2 ) 

plicity, we assume only one kind of outgas (water vapor) present, i .e., only R 

( 2 ) 

has non-zero values. The outgassing rote R (t) of crinkled sirtgle-oluminized 
mylar is shown in Fig. !• 

( 2 ) 

Ir our analysis, we use a set of R' '(t) values at different times successively. 

( 2 ) 

To obtain the solution at any chosen time, say t., we use the values of R' '(t .) 

0 ) 

and the chamber pressure p^' ^ (t^) in Eq$.( 19 )l Thus an appropriately chosen set 

of R^^^(t)ond p ^*^(t) will result in a solution of the P.(t), *the pressure history of 
c • 2 

the gos mixture inside the multilayer insulation system. However, the R (t) values 

are available only up to t = 48 hours; thus^our computations hove to stop there# 


I2 



Theoretical computation includes the pressures P. (t) inside each interstice, 
while only the pressure In the lost interstice Pj^(t) , i .e., the backside pressure, is 
experikiientally measured. Therefore, only the comparison for P^^ (t)can be made. 
Comparisons of the theoretical P^^ (t)and the experimental data are presented in 
Figures 2 through 6. 
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5« Conclusions 


First of all/ it should be noted that the comparison between the computed 
and measured values of the p^^'s presented in Figs. 2 through 6 rannoi Ih' taken 
seriously, since It is known tn.it tlie experiment for t lie dit erminat ion of tlie 
outgassing rate was not appropriately performed; tluis, tlie v.tliies of R''- (t) 
presented in Fig. I, upon wliieli our ealeulatlon was based, are not reliable.^ 
Wevertlieless, Figs. 2 tliroiigli 6 serve to siiow ;i qu.il i t.at i ve i-.omp.ir ison and 
to indicate some possible experimental errors. 


It can be seen from the figures that at early times the computed pj^ is 
much higher than the experimental values. However/ the computed and the mea- 
sured values are coming closer and closer together as time proceeds, 
and good agreement is established after a reasonably long time. Since tlie effect 
of preconditioning on the outgassing rate is likely to be a weak function of time 
when t is large/ one probable and important course of experimental discrepancy 

would arise from different preconditionings of samples/ i .e./ the samples For the 

( 2 ) 

measurement of (t)and the sample for the determination of ^(t) were dif- 
ferently pre-treated. Another source of error would be the contamination of 
equipments inside the vacuum chamber. 

It is believed that agreement between the theory ond experiment can be 
achieved if/ with special attention to preconditioning/ a set of consistent ex- 
periments is performed. 


At the time of the writing of this report/ the Marshall Space Flight Center is 
planning to rcconduct the outgassing rote experiment. 
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FIG* 1. Outgassing Rate for OinMed Single-Aluminized Mylar. 
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FIGt 3. Compaiiion of Computed and Mcasufod Backside Pressures/ 20-Laycr 
Sample . 


18 




in torrs 



FIG . 4. Comparison of Computed and Measured Baclcside Pressures/ 35"Layer 
Sample. 
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INTRODUCTION 


Heat exchangers are an essential part of the process of cooling various 
coHponents in space vehicles . Since* in any design related to space 
travel* imys of reducing ueight and volume are alvays of utmost ia^rtance* 
a concept that would lead to light and compact heat exchangers should cer- 
tainly be looked into. One possibility that has been suggested is the use 
of cryogenic liquids as the coolant In these heat exchangers. This offers 
an excellent possibility for improved performance and reduction in weight 
and volume since cryogenics provide extr»ely low t»perature& . This concept 
appears even more attractive when one considers the fact that cryogenics are 
already used as fuel in space vehicles and hence readily available to serve a 
dual purpose. However* when humid air comes in contact with a cryogenically 
cooled surface* frost forms rapidly on the surface and may eventually plug up 
the system and render it useless until the frost is removed. Therefore* the 
design of such a heat exchanger would require d<^tailed infoncatior. regarding 
this frost formation. This is where the purpose of this study comes in. 
Objectives ; 

The basic objective is to study the frost formation on super-cooled 
surfaces as a function of time* air velocity* air inlet and outlet tempera- 
tures* This would be done analytically as well as experimentally to arrive 
at empirical relations that would predict heat transfer coefficients and 
rate of frost growth for different conditions* 

Exp eriment al Se t-Up ; 

The study began with a crude experimental model (Figure 1) to gain a 
perspective of the problem* The set-up primarily consisted of a blower 
providing air flow for the heat exchanger* The heat exchanger consisted of 
two flat plates and these places were cooled by the cse of dry ice in an 



acetone bath« However, certain problems were realized immediately. First, 
the room air did not supply sufficient moisture to show a significant frost 
build up. Teats were conducted with Scott Phase Heat Exchanger (Hodel 9058) 
to investigate its use as a source of moisture for the system. This equip- 
ment proved to be successful in providing sufficient moisture to enable us to 
observe for the first time an actual plug-up of the air passage due to frost 
build-up. Photographs of this passage were taken at different intervals. 

This plug-up was observed for various width spacing of the two plates. In 
order to see how thick the frost muld get before It begins to act as an 
insulator, it would have been necessary to have a very large width between 
the plates and then maintain the humid flow and low temperatures for several 
hours. This was not possible for the existing set-up due to the problem of 
boll-off of the coolant and hence the inability to maintain cold plates for 
the required duration of time. 

Another problem that was realized was the fact that this set-up had no 
means of recording the frost growth. Thought was given to various methods 
of recording this frost growth. One such method consisted of measuring 
capacitance change between the plates and relating it to frost build-up. 
Hewlett Packard Universal Bridge Instruments (Model 4260A) was utilized for 
recording such capacitance changes. However, it was soon discovered that the 
capacitance between the plates changed not only due to frost build-up, but 
was also affected by the plexiglass casing, as w 11 as any changes in the 
flow rate or temperature of the air passing between the plates. It was 
therefore decided that this method was not feasible. 

Having realized that this set-up had no means of recording frost growth 
the air velocity, humidity and temperature, and that the existing test 
section allowed coo much coolant boil-off and had no provisions for draining 
the water from the melting of frost at the end of an experiment, it was 
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decided to design and assemble a new system and a new test ^.ection. A ( onsidi v- 
able amount of time was spent towards designing such a system, and the fin;'* 
arrangement is shown in Figure 2« 

Cov^ressed air is used as the source of flow for the system. This 
eliminates the need for a blower. A pressure regulator and a gate valve are 
used to control the flow rate. This flow rate is measured using a Cox Turbine 
meter and a frequency counter* Pressure and temperature are recorded at this 
point to convert from ACFM to SCFM. This air supply is now passed through an 
8" X 10" rectangular duct where it is heated to a desired temperature by a 
variable 5 kw duct heater (H.W. Tuttle Co.). Relative humidity is controlled by 
using a power spray humidifier and recorded on a Serdex Hygro thermo graph, 
before being sent to the test section. A picture of the assembled system is 
shown in Figure 3. 

The new test section has been designed but not yet constructed. A picture 
of the proposed test section is shown in Figure 4. When assembled, this test 
section will have the same two flat plates forming a rectangular passage. 
However, the new concept has an air passage below and above the test passage 
to circulate warm air eliminating fog formation and making it possible to 
record visually, froir above, the frost growth at various locations along the 
test section* Liquid Nitrogen will be used instead of dry ice to provide the 
cryogenic temperatures. (The boxes ci>iualn the spray-on foam insnlati<Mi 
(SOFI) to reduce boll-off.) Th is nt w lest sc‘cl ii>n h<is an easy a lor 

changing the spacing of Che two plates without dlaassembllng the whole thing. 

The handles that hold the plates in place travel in grooves and simply need 
to be moved In or out to change the width of the gap. 

Analytical Study : 

The analysis of an approximate model for Che growth of a depos*' ;ed layer 
on the cryogenic surface has been Initiated, the details of which have been 
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prevloualy reported In an Interim report Qtod, 1, Project II, October 
1972). 
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FIGURE 2. SCHEMATIC DIAGRA*-! OF A SYSTEM PROVIDING FLOW RATE, 
TEMPERATURE & HUMIDITY CONTROL OF AIR 
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FOREWORD 


The work described In this report is part of a continuing effor 
being carried out at UAH* The participants in this study have maintained 
a close liaison with Bob Fisher and Mr. Jerry Vanneman of NASA- 
MSFC, * The purpose of this close liaison has been to discuss 

problems t results » and new techniques in order that the results of this 
study be as beneticlaJ to this group as possible. Cbnsequentiy, many 
informal recommendations have been made» which may or may not be 
formalized in this report. 
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INTRODUCTION 


1.1 The Problem 

The problem to be studied stems from the requirement for t*'e control 
system thrustors of a space shuttle type vehicle, necessary for exo- .mospheric 
flight control, to be buried vithin the heat shield of the vehicle so that 
Interference heating will be reduced during flight in the atmosphere. The 
geometry of a typical installation is illustrated in Figure 1. During 
atmospheric flight, the non-firing thrustor nozzle is simply a cavity in 
the surface of the vehicle. The interaction of the hypersonic external flow 
%rith the cavity produces excess local heating vrhich must be ameliorated. 

Mechanical systems, such as a cover that slides over the nozzle exit, 
and active cooling systems introduce other problems besides adding undesirable 
weight and complexity to the vehicle. The simple expedient of increasing heat 
shield thickness in the region of the cavity would ultimately imply additional 
ablation in the area and a probable change in noz/. !c sii'.e or sh.ipo from flight 
to flight. This may require additional refurbishment before each flight. 

This study, therefore, is aimed at determining the mechanisms involved in 
producing the high local heating rates in and around the thurstor nozzle in 
non-firing condition. In such a case, the nozzle is treated simply as a cavity. 

1.2 Literature Survey 

A survey of the literature on this and related problems has revealed 
that very few results applicable tc this particular problem are available. Of 
the papers reviewed, only those referenced in this report appear to have direct 


consequence for this study. 



1.3 £xperln«ital Program 

A test program, conducted at AEDC by personnel of General Dynamics, Corp., 
has produced some results directly applicable to this study. In order to pro- 
vide a basij for comparison or for development of empirical relations, a short 
study vas made to compare the flow conditions in the wind-tunnel studies with 
flight conditions. Of particular concern was the question of whether the 
boundary layer ahead of the cavity was laminar or turbulent in either the 
wind tunnel or flight regimes. A map of these conditions is shown in Figures 

2 and 3. The wind tunnel conditions are for the von Karman Facility, Tunnel 
1 2 
B, at AEDC . The flight conditions are for a typical reentry corridor, with 

no special attempt to relate the reentry trajectory to the space shuttle 

vehicle. 

In the wind tunnel tests and in the flight regime, it is to be expected 

that the boundary layer ahead of the nozzle cavity will be fully developed 

and turbulent* In flight, but not in the wind tunnel, the boundary layer 

probably would contain significant amounts of ablation products. No attempt 

is made to account for the influence of ablation or surface heating ahead 

of the nozzle area, except to account for the total temperature profile through 

tne boundary layer and, hence, the total pressure profile. 

Since the analysis of the experimental data is being carried out at 

General Dynamics, no general discussion of the results will be made here. In 

so far as possible, comparisons of the results of this study are made with 

3 

the raw data of the experiments. These confirm the qualitative features 
of the flow derived from theoretical analysis. 
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ANAL\SIS 


2.1 Fundamental Features of the Flow 

In the flight vehicle, the thrustor nozzles are to be installed in 
curved surfaces. Presumably, the local radius of curvature of the body 
will be sufficiently large cf. the nozzle diameter so that the effects of 
surface curvature will be at most of secondary importance in this problem. 

The wind tunnel test program centered mainly on axlally**symr'etric nozzles 
set Into a flat-plate model. In this study, only one basic configuration 
is considered. This is the axlally-symmetric nozzle cavity set into a flat 
plate with the nozzle axis normal to the plate surface, as in Fig. 1. The 
flow would be two-dimensional in the absence of the cavity. The disturbance 
produced by the cavity Is three-dimensional, which greatly adds to the 
complexity of the problem. 

Sketches of a two-dimensional cavity flow and the central plane of the 
present problem are shown In Fig. A. In the two-dimensional problem, the 
boundary layer separates from the upstream edge of the cavity. The dividing 
streamline leaves the separation point and intersects the downstream surface 
of the cavity, at a stagnation point. The location of the stagnation point 
depends upon several parameters such as cavity height and cavity width. But, 
in the steady state 2-D solution, the stagnation streamline is the separation 
streamline. The flow within the cavity circulates about a vortex core, 
driven by the shearing stress along the dividing streamline. 

In the three-dimensional case, fluid may move in the y direction so that 
a net mass may be injected into the cavity in the central plane. In a steady 
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8tate» the mass flowing Into the cavity near the center line must flow 
out of the cavity at some distance from the central plane. The symmetry 
of the cavity clearly shows that stagnation should occur only in the 
central plane » so that the lateral flow would be expected to occur. The 
Injected mass flow would increa^^e the angular momentum of the flow in 
the cavity, thereby strengthening the vortex within the cavity. This 
vortex must exit the cavity with the outflow and be wat led downstream. 

The resultant vortex shape is reminiscent of the classical horseshoe 
vortex of a lifting wing of finite span. (Fig. 5) A crude model has 
been devised to indicate the influence of this vortex on the pressure 
distribution in the neighborhood of the cavity. The model is discussed 
In some detail in the next section of this report, ii eluding a comparison 
with experimental results. No attempt has been made to include com- 
pressibility effects, or the interaction of the trailing vortex with the 
boundary layer. 

The strength of the vortex within the cavity Is difficult to ascertain. 
Values may be inferred from two-dimensional results, but great care must be 
exercised in order to avoid erroneous conclusions. A discussion of two- 
dimensional and axially symmetric (annular cavity) studies is presented in 
a succeeding section of this report. As a result of these studies, some 
tentative conclusions may be drawn, particularly in reference to future 
experimentation* 

2.2 Horseshoe Vortex Model 

The vortex within the cavity in the three-dimensional case must extend 
”to infinity,” close on itself, or end on a solid wall, according to a theorem 
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due to Helmholtz. It has also been shown that in viscous flow» a vortex 
line may end on a solid boundary only at a stagnation point. Consequently, 
the vortex generated in the cavity in the case under consideration must 
have the form as described earlier (Fig. 5). 

The horseshoe vortex must be one of the dominant features of the cavlty-^ 
boundary layer Interaction. The flow Induced by the vortex would modify 
pressure distributions and velocity fields to a significant degree. Viscosity 
and compressibility would modify these effects, of course. 

A crude initial model was de^'ised by considering the superposition of 
a uniform flow on the flow due to a horseshoe vortex in steady, potential 
flow. Viscosity, compressibility and body forces are neglected. The 
boundary condition that there be no flow through the surface is Imposed by 
the method of Images. The cavity was not included in the geometry. 

The velocity field induced by a horseshoe vortex lying in an xy plane 

4 

at a height z above the surface is given by : 

i 



+ (x-x^)k 
(z-z,)^ + 




y -y. 




(z‘Zj) j + 6'}’y)<< 

+ (y ] y)^ 

(z-z j) j -t- (/j ” y)*^ 
(z-z^)^+ 


1 + 


X - X, 


^(x-x^)^ + (y-y^) +(z-z,) 



5 



where T is the circulation. The model is illustrated in Fig. 6. In com- 
pariBon of this initial model with the actual flow geometry, the length of 
segment A-B of the horseshoe vortex is taken as being equal to the nozzle 
exit diameter. 

Since the image system utilized above guarantees that there be no flow 
through the surface z • 0, this model shows no flow into the cavity. In 
order to obtain a dovnwash velocity over the cavity opening, an additional 
line vortex segment is superimposed, on the flow field in the upper half plane. 
This new vortex segment is parallel to the segment A£, at a distance equal 
to z, below the cavity surface. Thus the velocity vector is obtained: 


^ p 

\/ = 1 1 ; + — - 
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V U 1 ^ 
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In this model the compressibility, the boundary conditions on the cavity 
walls, and the boundary layer upstream of the cavity have not been taken 
into account* The pressure at any point is obtained from the incompressible 
Bernoulli Equation. 

2.3 Vortex Strength 

Since the velocity and pressure fields in the vortex model are depen- 
dent upon the strength of the vortex, F, (circulation) some method for 
predicting the value of F must be specified. In the absence of any theoreti- 
cal study which would specify F for the geometry of this problem, an appeal 
is made to the qualitative similarity between the flow in the plane of 
symmetry and the flow in a two-dimensional problem. 

Burggraf^ has given the following relation for the vorticity, o), in a 
two-dimensional rectangular cavity: 
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0.0615 (a/Xjj) 


2/31 1/2 


a |_(l+b/a) J 

where a is the width of the cavity, b is the depth of the cavity, is the 
"effective running length" to the separation point (in our case, the dis- 
tance from the leading edge of the plate to the edge of the cavity) , and 

— r 

F is the root-mean-square of the surface velocity over the perimeter of 
the cavity. 

The values of F* given by Burggraf are: 



2 

3 

.12488 

.15583 



.01092 .01732 .03122 I .07029 


For constant values of U,,,, Xq, and b, the variation of the vorticity w 
with cavity width, a, may be obtained from 


r — " i/2 

w /a\ (1 + b/Sj) F^* 

0 )j \a y (1 + b/a) F^ 


where the subscript 1 denotes a reference condition. With taken as 

being equal to b» this variation of u) with a is shown in Fig. 8. 

The heat transfer rate on the reattachtrent (downstream) wall of the 
cavity is given by Hodgson^ as: 


This variation is also shown in Fig. 8. 

Burggraf *8 theory is based upon the plane, Invlscid, incompressible, 
but rotational, flow over a rectangular cavity. According to his theory. 
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the heat transfer rate depends only on the cavity width-to-depth ratio. 
Hodgson includes compressibility. His results indicate that the heat 
transfer rate may vary with cavity width, even with b/a fixed, thus re- 
vealing a scaling effect. 

A comparison of the heat transfer rates on the reattachment wall 
predicted by these two-dimensional theories with the rates obtained in the 
plane of syimivetry in the three-dimensional tests in a wind tunnel are 
shown in Fig. 9. Two items are immediately Indicated by this comparison: 

1) the test data is much steeper in the reattachment zone than that pre- 
dicted by either Burggraf*s or Hodgson's theory; and 2) the tests show a 
decrease in reattachment heating with decreasing nozzle diameter as evi- 
denced by tests with a two inch diameter nozzle, end tests with a one Inch 
diameter nozzle. Hence, scaling may be a very important factor in the appli- 
cation of the wind tunnel test results to the flight vehicle. 

The vortlclty predicted by Burggraf's two-dimensional theory was taken 
to predict the circulation F of the horseshoe vortex in the three-dimensional 
model discussed earlier. The pressure distribution obtained from this 
model are shown in Fig. 10. Experimental results, shown in Figs. 11-21, 
show qualitative agreement. 

2.4 Experimental Results 

Tests of the axially symmetric nozzle cavity in a flat plate at various 
angles of attack, were conducted in the AEDC-VKF Tunnel B, for various 
Reynolds Nurbers and at a nominal Mach Number of 8. Wind tunnel operating 
conditions for the various tests are listed in Table I. The flow conditions 
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on the flat plate for the ten cases shown in this report are listed in 
Table II. Isobars and isotherms in the vicinity of the cavity for these 
ten cases are shown in Figs. 11-21. 

The major qualitative features of the pressure field in these tests 
are predicted by the horseshoe vortex model of the flow. Of particular 
interest in this respect is the high pressure zones just outside the cavity 
at 0 “ 0^. Where the theory shows a decreasing pressure, the experiment 
shows a low pressure on the cavity lip and a rising pressure downstream. 

This discrepancy could be due to separation as the reattached flow tries 
to negotiate the sharp lip of the cavity exit. In any event, the agreement 
is remarkable in view of the crudeness of the theoretical model. This 
clearly demonstrates that the horseshoe vortex la the key feature of the 
flow field. 

One aspect of the test results Is very Interesting. In an effort to 
reduce the peak heating rate In the cavity, which occurs at the reattach- 
ment point, cold gas was bled (l.e., at low thrustor chamber pressure) 
through the nozzle in cases 2 and 3. With increasing bleed rate the stag- 
nation point (reattachment point) moved farther into the cavity and the 
reattachment wall temperature increased. Thus the nozzle bleed aggravated 
the heating problem instead of ameliorating it. The horseshoe vortex model 
of the flow explains this phenomenon. Since for the bleed cases (as opposed 
to thrustor firing) the stagnation pressure of the bleed gas is much, much 
less than the stagnation pressure of the external flow. As indicated In Fig. 
22, the recirculating gas Is not much Influenced by the bleed gas which 
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must separate from the downstream (reattachment) wall just above the nozzle 


throat. Consequently, all of the bleed gas must be directed upward along 
the rostream wall. This Injection of the bleed gas i ncreases the angular 
momentum of the cavity flow> and therefore Increases the vortex strength. 

The stronger vortex, in turn pulls in more of the external flow, again 
increasing the angular momentum in the cavity. The new stagnation stream- 
line, moreover, originates higher in the upstream boundary layer so that it 
has a higher stagnation temperature than the streamlines closer to the wall, 
due to heat transfer to the plate. 

2.5 Cavity Shape 

Thomke^ has studied the flow of turbulent boundary layers flowing over 

axially-sytnmetrlc cavities (annular grooves cut into the surface of a right- 

circular cone). The experimental tests were conducted at Mach Ntmbers of 

2.0, 2.5, 3.0, and 4.0 and wind tunnel unit Reynold’s Numbers of 0.95, 1.02, 

1.36, 1.43, and 1.30 million. A comparison of pressure distributions taken 

atM“2*0, R «0.95x 10^ per inch for two cavity shapes is shown in Fig. 
e 

23. 

The test results clearly show that rounding the reattachmert wall pro- 
duces a much flatter (almost constant) pressure distribution along the cavity 
floor. Additionally, the rounded corner results ti a lower peak pressure, 
and the point of maximum pressure being moved downstream. These results 
would imply that rounding the corner would reduce heating in or near the 
cavity. 
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APPLICATION OF 2-D THEORY 


3.1 Modified 2-D Theory 

As Indicated before, the flow over the actual thrustor nozzles and over 
the axlally-symmetric nozzles mounted in the wind tunnel models is three- 
dimensional. But presently, there is no theoretical stud> on the three 
dimensional cavity flow which can be used to predict the heat transfer rate 
to the thrustor nozzle walls, fince, due to symmetry of the cavity, the flow 
would stagnate in the central plane, the central plane would have the most 
severe heating caused by the flow. Further, the flow in the plane of symmetry 
is qualitatively similar to a two-dimensional flow. Hence, an attempt was 
made to utilize the existing two-dimensional theory coupled with the 
experimental data, to predict the reattachment heating on the central plane 
of the thrustor ./zzles. 

From the literature survey, it was found that Hodgson’s^ theory was the 
most applicable to this study* The model chosen in Hodgson's theory is shown 

o 

in Fig. 24. It is similar to that used by Chung and Vxegas,^ who neglected 
the oncoming boundary layer. The resulting expression for the velocity down 
the reattachment wall is 

• 1 - 0.189 fexp (-5.3 f ) 

“d 

_ 10.6 I exp [-y/L (9.87 + 28.1) ^ ] 

n*! 2 

9.87 n + 28.1 (6) 
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The effect of the upstream boundary layer is accounted for through Its 

influence on the dividing streamline velocity U , which can be obtained 

d 

from the analysis of Denison and Baum.^ A least-square curve fit of 
Denison and Baum’s result has been given by C. D. Engel as 

^d 

- a + a C + a C* + a + a C** + a r® (7) 

Wg 012 3 4 5 

where 

S* » w/x (flat place) 

0 

w ■ cavity width 

X ■ effective* running length 
c 

c •= logjo (S*) 

a - 0.: 1706 
0 

a » 0.088707 
1 

a —0.049822 
2 

a —0.0020116 

3 

a - 0.0047757 

4 

a - 0.0006549 

5 

The expression for the local Stanton number is 



for air (p^ » 0.7), Cj - 0.413, c^ - 0.435, c^ - 1.87 and the dividing 
streamline temperature is found from Schllchtlng.^^ 
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The cavity temperature vat. assumed as 


T - 0.5 (T + T ) 

c d w 


( 10 ) 


eliminate T from Eq. (9) and (10), expression of T is 

d 



( 11 ) 
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The reattachment zone length L was given by Neater, 


7-1 


L ■ l'^ 


1 + .447 2 


5.22 + 4.41 - 1 

\T^ i 


where is the value corresponding to incompressible flow. 
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with above formulas, the heat transfer rate h is giver by: 


St V c p (T , - T ) 

y paw 

- T ) 

0 W 


(14) 
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The above nodal and f^rvulaa (except Eq. (11)) were obtained from Hodgson's 
theory for a tiio-dlneas tonal rectangular cavity. In order to apply this 
theory to predict the heat transfer rate in the central plane of axiallv- 
aynnietric nozale» simple modifications based upon empirical data were 
made, leading to the folioving formulae: 



(15) 



(16) 


Equation (15) accounts for the flow reattaching inside the cavity instead 
of reattaching right at the comer, as in 2-D theory. Equation (16) 
describes the heat transfer to the cavity vail, strongly influenced by 
the horseshoe vortex within the cavity. 

3.2 Zero Incidence Cases 

A calculation, based on the above formula, wa<. ^de to compare with 
AEDC test data as shown in Fig. 25. For zero incidence cases, Chi<^ 
modified Hodgson's method matches AEDC test data pretty well. This 
strongly indicates that the flow before the cavity was laminar during 
the test of zero incidence cases. This could happen because during the 
test of zero incidence cases, a slightly negative angle of attack was made 
".o compensate for the displacement effect of the boundary layer. And 
sicr.e the growth of boundary layer is non-linear, by doing tnis com;jensation 
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the boundary layer vould tend to be stabilized because of a slightly 
favorable pressure gradient a%iay from the leading edge. Hence, in spite 
of the high Reynold's nuaber in the test condition, the flow icay well be 
stabilized sufficiently that the onset of transition occurs downstream of 
the cavity. 

3.3 Angle of Attack Cases 

In Fig. 27, in the angle of attack cases, there is an order of magni- 
tude increase in the normalized experimental heat transfer rate above the 
zero angle of attack results. This increased heating is not predicted by 
the model described above* In order to find the reason for this order of 
magnitude discrepancy, several lices of study were followed: 

(a) External Flow Properties - The pressure data shoved that there 
was a weak shock ahead of the first point of the test data, which 
was 7«65" from the center of the cavities. This shock might be due 
to the Junction of the test plate with the leading edge* This 
Junction was 8.76” from the center of cavity. But calculation showed 
that this weak shock did net have much influence on the external 
flow in the vicinity of the cavity and would nor thereby effect 
reattachment wall heating rate. 

(b) The plot of actual heat transfer rate. Fig. 28, showed that 
the order of magnitude discrepancy was in the actual heat transfer 
rate and not due to a poor choice of h 

ref 

(c) Since the parameter which characterized the Hodgson's theory js 
the velocity along the dividing streamline, a calculation, which did 
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not include the empirical factor of 0*5 in £q. (16), based on the 
assumption of laminar flow ahead of the cavity but with increased 
dividing streamline velocity, which corresponding the maximum heat 
transfer rate, was carried out. The results. Fig. 29, are comparable 
to the test data* 

(d) Since in the test data there was no clear Indication of whether 
the flow before the cavity was laminar or turbulent, a calculation 
of the heat transfer rate in the central plane based on the assumption 
of turbulent boundary layer ahead of the cavity was also carried out. 
The dividing streamline velocity in the turbulent shear layer was 
given by Lamb.^^ 

The results of this calcilation are shown in Fig. 30. These results 
were about the same as that of the maximum rate possible in laminar 
flow as calculated in part (c) and were comparable to the test data. 

It should be noted that in this calculation only the dividing stream- 
line velocity was calculated for a turbulent shear layer, whereas 
both the velocity in recompression zone (See Fig. 24) and the local 
Stanton number were still calculated by using laminar flow equations. 
It ia aiso worth noting that only a change of laminar boundary 
layer to turbulent boundary layer before the cavity has brought the 
heat transfer rate to the maximum possible value in laminar flow. 

This reveils the importance of the study of turbulent flow. Hence, 
a more thorough study based on the consistent turbulent flow is needed 
for future work. 
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RESULTS 


Several corxlusions may be drawn from the preceding discussion: 

1. The boundary layer-cavity Interaction generates a 
horseshoe vortex reminiscent of the vortex system 
of a wing of finite span. 

2. The pressure distribution in and around the cavity 

is influenced very strongly by this horseshoe vortex. 
The influence of this vortex will extend far down- 
stream, and laterally beyond the sides of the cavity. 

3. Bleeding cold gas through the nozzle increases the 
strength of the vortex and, consequently, increases 
the peak heating rate in the nozzle. 

4. The key to passive system heat protection in the 
region of the cavity is the control of the vortex 
strength by shaping of the nozzle. Variation of 
cavity geometry could tend to weaken the vortex, 
thereby relieving the heating problem. 

5. For zero incidence cases, the method modified from 
2-D laminar theory gives good prediction in the 
central plane of axial ly-symmetric nozzle. 

6. For angle of attack cases, based on the assumption 
of turbulent boundary layer ahead of the cavity, the 
prediction of heat transfer rate is comparable to 
AEDC test data. 
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7. The heat trar*)fer rate on the reattachment wall 

calculated based on the turbulent velocity profile 
was an order of magnitude higher than that calculated 
based on laminar flow. 
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Table I. Wind Tunnel Operating Conditions From 
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Table 11. Fxac Plate Flow Field Conditions in A£DC Test 
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rig* 2. Comparison ol Mach Numbers between Re-entry Flighi and 
Wind Tunnel Test Conditions 


25 




Re/ ft 



and Wind Tunnel Teat Conditions 


26 







a) 



b) 


FIG . 4. Sketches of Flow in Cavity • 

a) Two-dimensional Flow 

b) Rhine of Symmetry in 3-D Cose 
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FIG - 6- Superposition of Horseshoe Vortex ond Free Streom, Using Method of 
Images to Enforce no Flow Through the Surfoce- 
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FIG 7 . Horseshoe Vortex Model With Additioral Vortex Segment to Induce 
Flow Through Cavity Fxit Plane. 
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TEMPERATURE ( lower half, ®R ) DISTRIBUTIONS 
INSIDE AND AROUND THE EDGE OF CAVITY, UNDER 
CASE 1 




35 


zed by P-inf ) AND 











36 


INSIDE AND AROUND THE EDGE OP CAVITY 
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Flg« 16 . PRESSURE ( upper half, nonnallzed by P-inf ) AMD 

temperature ( lower half, **R ) DIS1RIBUTI(3NS 
INSIDE AMD AROUND THE EDGE OF CAVITY, UNDER 
CASE 3 
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Fig* 18. TEMPERATURE DISIRIBUTIOH (*R) IHSIDE AND ARCWND 
TRE edge of the CAVi."T, UNDER CASE 4 
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Fig. 20. TEMPERATURE DISTRIBUTIGH (®R) INSIDE AND AROUND 
riE EDGE OF THE CAVITT, UNMH CASE 5 
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Flg.2l. PRiSSSURE DISIRIBaTION ( upper half, normallxed 
by P-inf ) AND TEMPERATURE DISTRIBUTION ‘*R ( 
lower half ) UNDER CASE 6 
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Fig« 2S« Comparison between modified Hodgson's thei and test 
data for zero Incidence cases 





^ te$t data from Nlcoll 
R^/ln =5.9*10^ 
cone angle 2(f 
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value from modified Hodgson* s theory 
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Heat Transfer Bate q (Btu/ft-sec) 






Pig. 29 Comparison between Modified Hodgaon's Theory (nith increased 
di'^idlng streamline velocity) and test data foi: angle of 
attack cases 
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